Defect energy of infinite-component vector spin glasses 
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We compute numerically the zero temperature defect energy, AE, of the vector spin glass in the 
limit of an infinite number of spin components m, for a range of dimensions 2 < d < 5. Fitting 
to AE ~ L° , where L is the system size, we obtain: 9 ~ — 1.54 = 2), ~ — 1.04(d = 3), 
8 ~ —0.67 (d = 4) and 8 ~ —0.37 (d = 5). These results show that the lower critical dimension, di, 
(the dimension where 8 changes sign) is significantly higher for m = oo than for finite m (where 
2<di< 3). 
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I. INTRODUCTION 

There has recentlji 2 ^ been interest in spin glasses 
where the number of spin components, m, is infinite, be- 
cause this limit provides some simplifications compared 
with Ising (m =1), XY (m = 2), or Heisenberg (m = 3) 
models. For example, in mean field theory (i.e. for the in- 
finite range model) there is no "replica symmetry break- 
ing"— so the ordered state is characterized by a single 
order parameter q, rather than by an infinite number 
of order parameters (encapsulated in a function q{x)) 
which are needed^ for finite-m. In addition, there are 
special numerical techniquesi^iiSii which can be used 
to study finite-range m = oo spin glasses in which the 
(finite) sample is solved exactly without the statistical 
errors and equilibration problems inherent in the Monte 
Carlo methods used for finite-m. 

There are, however, significant differences between 
Ising, XY, and Heisenberg spin glasses, on the one hand, 
and m = oo spin glasses on the other. In the Ising spin 
glass in three dimensions there is clearly a finite tem- 
perature transition 8 ., and we have argued* that the same 
is true for XY and Heisenberg spins, though the latter 
is still somewhat controversial, see e.g. Refs. TiTMTTirT^ . 
Hence, for m = 1,2 and 3 the lower critical dimension, 
below which T c is zero, is less than 3 (in fact 2 < di < 3). 
However, for m = oo, one finds^i T c = in three- 
dimensions, so di must be greater than 3 in this case. In 
fact, VianaAi make the surprising claim that di = 8 for 
m == oo, by attempting to sum up the perturbation ex- 
pansion. Curiously the upper critical dimension (above 
which the critical exponents have mean field values) is 
also predicted^ to be d u — 8, which is again different 
from the value for finite m where d u — 6. 

In this paper we attempt to determine the lower crit- 
ical dimension of the m = oo spin glass by computing 
the zero temperature "defect energy" , AE, for a range of 
dimensions, 2 < d < 5. The defect energy is the charac- 
teristic energy change when the boundary conditions are 
changed from periodic (say) to antiperiodioi£*i£*iZi. It is 
expected that 

AE ~ L e (1) 
where L is the system size and 9 is a "stiffness exponent" . 



If > then the system is stiff on large lengthscales and 
so one expects T c > whereas if 9 < then it costs very 
little energy to break up the ground state configuration 
at large scales, and so presumably T c = 0. Hence di is 
the dimension where 9 — 0. For the case where 9 < 0, so 
T c = 0, the correlation length £ diverges as T — > like 
£ T~ v , and standard scaling argumentsi^i^ then show 
that v = -l/d. 

Our main result is that 9 < for the full range of 
dimensions (2 < d < 5) that we are able to study showing 
that di is significantly greater than 5, i.e. much larger 
than for finite-m. 

In Sec. II we discuss the models and numerical im- 
plementation. In Sec. Ill we discuss our results from 
simulation in two to five dimensions. We give our con- 
clusions in Sec. IV. 



II. MODEL AND METHOD 



We take the Edwards- Anderson 18 Hamiltonian 



1~t — y JijSi • Sj , 
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where the spins Si (i = l,---,N) are classical vectors 
with m components and normalized to length m 1 / 2 , i.e. 
S 2 = m. The summation is over nearest-neigbor pairs. 
The interactions Jij connect nearest neighbors and are 
independent random variables with a gaussian distribu- 
tion with zero mean and standard deviation unity. 

At finite temperature and for m = oo the spin-spin 
correlation functions, 
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are obtained from the following set of equations 
T~ 1 C l j = (A- 1 ).., where 



1,2,3,6,7 



(4) 
(5) 



Here the Hi (i = 1, • • • , N = L d ) are Lagrange multipliers 
enforcing the normalization of the spins: 
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(i = l,2,---,N). 
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To proceed, one solves the N equations, Eq. JfjJ), to obtain 
the Hi, and then determines the Cy from Eqs. (@J and 

©■ 

At zero temperature, Eqs. J1J and Eq. 101 are no longer 
well defined. However, since there are no thermal fluctu- 
ations, each spin lies parallel to its local field: 



(7) 



Remarkably, it was shown by Hastings^ that these local 
fields are precisely the zero temperature limit of the Hi 
in Eq. (J5J. Another interesting result found by Hast- 
ings is that the number of independent spin compo- 
nents mo used to form the groundstate satisfies a bound 
mo < y/2N. By independent we meant that we can al- 
ways define coordinates for the spins such that the pro- 
jections of the spins are only non-zero for mo directions 
and no spin components are found in remaining to — too 
directions. Furthermore, it is found numerically that 



N 1 *, (to < to), 



(8) 



where [• • -] av denotes an average over disorder, and the 
values of /i we find for dimensions between 2 and 5 are 
given m Tabled 

We now see that to study the to — > oo limit, we simply 
need that to should be greater than mo- Since calculating 
Too involves finding the eigenvalues of a large matrix (N x 
N) for N = L d , it is computationally intensive for large 
L and d. In practice, we can determine [mo] a v from a 
smaller range of sizes and fit to Eq. ©. This allows us 
to extrapolate the value of [TOo] av to larger sizes. We then 
choose to to be significantly greater than [m,o] a v for all 
L. Note some tolerance is required because the precise 
value of Too varies from sample to sample. 

To find the groundstate, we use a "spin-quench" 
method^. Firstly, the local field on Si is computed via 
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,1/2 



(9) 



Next we set Si according to Eq. (0). This procedure 
is applied to each spin of the lattice sequentially until 
convergence, and then iterated to convergence. Our con- 
vergence criterion is that the magnitude of the change 
in each spin is less than about 10~ 7 . For finite to, this 
method does not guarantee the groundstate as there are 
many solutions to Eq. Q. However, it works in the 
to — ► oo limit because there is a unique solutionis in 
that case. 

To determine the defect energy AE, we first find the 
groundstate energy with periodic boundary conditions, 
E p , for a given set of bonds. Next, we reverse the L d ~ x 
bonds which wrap around the system in one direction (x 
say), i.e. the bonds that connect sites with x = 1 to 
those with x — L. We then obtain the groundstate en- 
ergy, E a , for these "anti-periodic" boundary conditions. 
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0.29 


-1.54 ±0.02 


0.65 ±0.01 


0.65 ± 0.05 


0.65 ± 0.02 
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0.33 


-1.04 ±0.02 


0.96 ± 0.02 


1.23 ±0.13 


1.01 ± 0.02 
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0.35 


-0.67 ±0.04 


1.49 ± 0.09 




1.5 ±0.1 
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0.37 


-0.37 ±0.07 


2.70 ±0.51 







TABLE I: Estimates of n and 9 — —1/v. The values of v are 
compared to estimates from finite temperature simulations 3 
and a previous calculation 7 at T = 0. 
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TABLE II: Number of samples used in the defect energy 
calculations. 



On average, neither periodic or anti-periodic boundary 
conditions is preferred and so we average the absolute 
value of E p — E a over many different configurations of 
bonds, i.e. the defect energy is defined to be 

AE= [\E P -E a \] m . (10) 

We expect that AE scales with L according to Eq. (JTJ. 



III. RESULTS 

We have performed simulations for dimensions d =2, 
3, 4 and 5. The number of samples for each size and 
dimension is presented in Table [H] The defect energies 
for different dimensions are plotted in Figs I1I2R1 and 01 
together with the fit to Eq. QJ. To reduce finite size 
effects when doing the fits, for d = 2 we omit the two 
smallest sizes (L = 4 and 10), and for higher d we just 
omit the smallest size (L = 4). 

The exponent 9 thus obtained is shown in Tabic [I] 
together with that obtained from previous calculations. 
The results for d = 2, 3 and 4 agree well with those of 
Ref. though we have better statistics than in that work 
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FIG. 1: d = 2: graph of AE against L with 9 = -1.54. 
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FIG. 3: d = 4: graph of AS against L with 9 = -0.67. 
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FIG. 2: d = 3: graph of AS against L with 9 = -1.04. 
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FIG. 4: d = 5: graph of AE against L with 9 = -0.37. 



and cover a larger range of sizes (128, 24 and 12 as op- 
posed to 12, 7 and 5). We are not aware of any other 
results for d = 5. Comparing with Ref. |sj, our results for 
d = 2 agree very well, while those for d = 3 are a little 
different, at the level of about 2a, which may reflect some 
corrections to scaling. 

In Fig. we plot our results for 9 as a function of 
d, together with a smooth curve through the points. It 



is obviously desirable to know the dimension, d\, where 
(9 = 0, but extrapolation of our data to larger d is very 
uncertain. However, it is clear that di must be signifi- 
cantly greater than 5, and hence much greater than the 
its value for finite m, which is between 2 and 3. It is not 
possible to test precisely the claim of Viana^ that di = 8, 
because we can't estimate 9 for d close to 8. However, 
our data does not rule out this possibility. 
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FIG. 5: Graph of 6 against d. The solid line is a smooth 
curve through the data. Its extrapolation to d > 5 (curve 
with long dashes) is, however, very uncertain. 



IV. CONCLUSIONS 

We have computed the zero-temperature stiffness ex- 
ponent 9 for the vector spin glass in the limit where the 
spins have an infinite number of components. We have 
obtained better statistics for a larger range of sizes and 
dimensions that in previous work£. Our results for 9 
agree with those of Ref. for the dimensions (d = 2- 
4) considered by them, while our result for d = 5 is new. 
The trend in our data, shown in Fig. indicates that the 
lower critical dimension must be significantly larger than 
5, the largest dimension we have been able to study, and 
may equal 8, as predicted by Viana^, but it is currently 
not technically possible to determine 9 for sufficiently 
high dimension to test this claim precisely. 
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